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Abstract We investigate the excitation of magnetoacoustic-gravity waves gen- 
erated from localized pulses in the gas pressure as well as in vertical component 
of velocity. These pulses are initially launched at the top of the solar photosphere 
that is permeated by a weak magnetic field. We investigate three different config- 
urations of the background magnetic field lines: horizontal, vertical and oblique 
O [ to the gravitational force. We numerically model magnetoacoustic-gravity waves 

by implementing a realistic (VAL-C) model of solar temperature. We solve two- 
dimensional ideal magnetohydrodynamic equations numerically with the use of 
the FLASH code to simulate the dynamics of the lower solar atmosphere. The 
initial pulses result in shocks at higher altitudes. Our numerical simulations 
reveal that a small-amplitude initial pulse can produce magnetoacoustic-gravity 
waves, which are later reflected from the transition region due to the large tem- 
perature gradient. The atmospheric cavities in the lower solar atmosphere are 
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found to be the ideal places that may act as a resonator for various oscillations, 
including their trapping and leakage into the higher atmosphere. Our numerical 
simulations successfully model the excitation of such wave modes, their reflection 
and trapping, as well as the associated plasma dynamics. 

Keywords: Waves, Magnetic fields, Corona, Granulation 
1. Introduction 

The complicated magnetic field configuration of the Sun plays a key role in 
various types of dynamical plasma processes in its atmosphere, including all 
the significant plasma dynamics of the lower solar atmosphere. The resulting 
magnetic structures channel energy from the photosphere into the upper at- 
mosphere, in form of magnetohydrodynamic (MHD) waves, and such waves 
experience mode conversion, resonances, trapping and reflection which results in 
the complicated dynamical processes in the lower solar atmosphere, the details of 
which depend on the plasma properties as well as strength of the magnetic field. 
The complex magnetic field and plasma structuring in the lower solar atmosphere 
support the excitation of various kinds of MHD waves, and their propagation, 
reflection and trapping has been studied extensively both on theoretical and 
observational grounds (e.g. Murawski et al., 2011; Srivastava & Dwivedi, 2010; 
Srivastava, 2010; Fedun et al., 2009; Srivastava et al., 2008; Hasan et al., 2005; 
McAteer et al., 2003; Gruszecki et al., 2011, and references therein). The evolving 
magnetic fields of the lower solar atmosphere also lead to transient processes 
across a wide range of spatial-temporal scales in form of the eruption and 
associated phenomena. For example, various types of solar jets are formed at 
short spatial-temporal scales which play a significant role in mass and energy 
transport and also couple the various layers of the solar atmosphere (e.g. Shibata 
et al., 2007; Katsukawa et al., 2007; Srivastava & Murawski, 2011, and references 
therein). In addition, the magnetic activity and injections of helicity into the 
lower solar atmosphere result in large-scale eruptive phenomena, including solar 
flares and coronal mass ejections (CMEs) in the outer part of the magnetized 
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solar atmosphere (e.g. Srivastava et al., 2010; Shibata & Magara, 2011; Zhang 
et al., 2012, and references therein). Therefore, the coupling of the complex 
magnetic field in various layers of the Sun due to waves and transients is one of 
the most significant areas of contemporary solar research. 

In the quiet-Sun magnetic networks, cavities are important locations where 
the magnetic fields are sufficiently inclined due to their well evolved horizontal 
components. They are formed over the granular cells in form of the field-free 
regions due to the transport of plasma at their boundaries and overlaid by 
bipolar magnetic canopies. The vertical magnetic fields, however, reside mostly 
in the core of such magnetic networks (Schrijver & Title, 2003; Centeno et al., 
2007). Such magnetic cavity-canopy systems are thought to be ideal resonators 
of the various MHD waves that can be trapped in the cavity, as well can also 
leak in form of the magnetoacoustic-gravity waves upward through the core of 
such magnetic networks. It is thought that the field- free cavity regions under- 
lying the bipolar canopy can trap the high-frequency acoustic oscillations, and 
the low-frequency components may leak into the higher atmosphere in form 
of magnetoacoustic-gravity waves (Kuridze et al., 2008; Srivastava et al., 2008; 
Srivastava, 2010). Therefore, such magnetic structures in the lower solar atmo- 
sphere are the regions that may play an important role in the wave filtering (e.g. 
Mcintosh & Judge, 2001; Krijger et al., 2001; McAteer et al., 2002; Vecchio et 
al., 2007, 2009; Srivastava, 2010, and references therein). 

Alongside recent high-resolution observations of MHD waves in the lower 
solar atmosphere, extensive efforts have been made in the area of analytical 
and numerical modeling of such waves: Fedun et al. (2009) have investigated 
the 3D numerical modeling of the coupled slow and fast magnetoacoustic wave 
propagation in the lower solar atmosphere; and recently Fedun et al. (2011a) 
have reported the first numerical results of the frequency filtering of torsional 
Alfven waves in the chromosphere. Apart from general numerical modeling of the 
waves in the lower solar atmosphere, models have also investigated the acoustic 
wave spectrum in the localized magnetic structures of the lower solar atmosphere, 
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e.g., magnetic cavity-canopy system (e.g. Kuridze et al., 2008; Srivastava et al., 
2008; Kuridze et al., 2009, and references cited therein). 

It is also noteworthy that Bogdan et al. (2003), Fedun et al. (2009) and Fcdun 
et al. (2011a) discussed in detail the excitation, propagation and conversion of 
magnctoacoustic waves in a realistic 3D MHD simulation. However, in these 
references, the waves driven by a periodic driver were discussed, whereas we 
numerically simulate the excitation of magnctoacoustic-gravity waves generated 
by pulses in the gas pressure and vertical component of velocity, mimicking an 
isolated solar granule. We aim to investigate and understand this simpler (albeit 
complex enough) system before tackling the more realistic, multiple granule 
system. Our philosophy is to build up our models incrementally, with a clear 
focus on the underlying physical processes at each step. 

In this paper, we investigate the excitation of magnetoacoustic-gravity waves 
generated from localized pulses in the gas pressure as well as in vertical compo- 
nent of velocity, modelling the effect of an isolated solar granule. These pulses are 
initially launched at the top of the solar photosphere that is permeated by a weak 
magnetic field. We investigate three different configurations of the background 
magnetic field lines: vertical, horizontal, and oblique to the gravitational force. 
We aim to show that small amplitude perturbations that are associated with 
such a granule are able to trigger large amplitude, complicated oscillations in 
the solar corona, which exhibit periodicities within the detected range of 3 — 5 
min. 

The structure of the paper is as follows: In Sect. [5] we describe the numerical 
model. We report the numerical results in Sect. [3] and present the discussion and 
conclusions in Sect. [4] . 



2. Numerical model 

We consider a gravitationally-stratificd solar atmosphere that is described by 
the ideal two-dimensional (2D) MHD equations: 

§+V-(eV) = 0, (1) 



SOLA: KM-AKS-JAM-R0_Solar_Phys_ver2e.tex; 30 August 2012; 0:11; p. 4 



Magnetoacoustic-Gravity Waves in the Solar Atmosphere 



dV 1 
g- + ? (V.V)V = -V P +-(VxB)xB + ffi , (2) 

£| + V-( P V) = (l- 7 )pVV, (3) 

-=Vx(VxB), V-B = fl. (4) 

Here g is mass density, V is the flow velocity, B is the magnetic field, p = 
k^gT/m is the gas pressure, T is the temperature, 7 = 5/3 is the adiabatic 
index, g = (0, — g, 0) is the gravitational acceleration, where g = 274 m s -2 , 
m is the mean particle mass and fee is Boltzmann's constant. Throughout 
this paper, we use the Cartesian coordinate system with the vertical 
axis denoted by y and the horizontal axis x. Henceforth, we assume 
that the medium is invariant along the z-direction with d/dz = and 
set the ^-components of plasma velocity and magnetic field equal to 
zero, i.e. V z = and B z = 0. The latter assumption removes Alfven 
waves from the system but still allows magnetoacoustic-gravity waves 
to propagate freely. 

2.1. Equilibrium configuration 

We assume that at its equilibrium the solar atmosphere is static (V e = 0) and 
threaded by a straight magnetic field, 

B = B §, (5) 

where s is a unit vector that is either vertical, horizontal or oblique. We choose 
-Bo by requiring that at the reference point (0, 10) Mm, the Alfven speed 

c A (y) = -j== (6) 

and sound speed 



Cs(y) lw (7) 

satisfy the following constraint: c\(y = 10 Mm) = 10 c s (y = 10 Mm). This 
constraint reproduces typical conditions in the solar corona where, typically, 
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c s = 0.1 Mm s 1 and ca = 1 Mm s 1 . As a result, the solar corona is magnetically 
dominated with B — 14.5 x 10~ 4 T. The plasma /?, 



M \ 7 c 2 s {y) 



(8) 



in the solar corona attains a value of (3(y = 10 Mm) = 0.012. It grows slowly with 
depth within the chromosphere and below y = 2 Mm reaches abruptly a value 
of (3{y = OMm) ~ 4 x 10 5 at the bottom of the photosphere (Fig. [TJ bottom). 
This large value of /3 evidences that in these low regions of the solar atmosphere 
the effect of magnetic field is negligibly small. 
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Figure 1. Equilibrium profile of solar temperature (top-left), mass density (top-right), 
and the plasma /3 (bottom). 

The straight magnetic field of Eq. (J5J) is a simplified model of curved magnetic 
field lines which are located in low regions of the solar atmosphere. However, as 
the plasma (3 grows fast with depth below y = 2.5 Mm reaching large values 
there (see Fig. [TJ bottom) we actually model the very weakly-magnetized low 
layers of the atmosphere, which do not correspond to flux-tubes and sunspots. 
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Therefore, in this first approximation, the curved magnetic field lines can be 
replaced by a straight magnetic field, and our models are justified. 

In Eqs. (|6]) and ([7|l. g e {y) and p e {y) denote equilibrium mass density and gas 
pressure, respectively. They are specified by the hydrostatic constraint which in 
the context of Eq. (JSJ states that the pressure gradient is balanced by the gravity 
force, 

- V Pc + g c g = 0. (9) 
With the ideal gas law and the y-component of Eq. © , we arrive at 



Pe(y) = Po exp 

where 



y dy 



A(y') 



, \ Pe(y) , 1fV . 



m - (id 

mg 

is the pressure scale-height, and po denotes the gas pressure at the reference level 
that we choose in the solar corona at y r = 10 Mm. 

We adopt an equilibrium temperature profile T e (z) for the solar atmosphere 
that is close to the VAL-C atmospheric model of Vernazza et al. (1981): see 
Fig. [TJ left-top panel. Note that T c attains a value of about 5700 K at the top 
of the photosphere which corresponds to y = 0.5 Mm. At higher altitudes T c (y) 
falls off until it reaches its minimum of 4350 K at the altitude y ~ 0.95 Mm. 
Higher up, T e (y) grows gradually with height up to the transition region which 
is located at y ~ 2.7 Mm. Here T e (y) experiences a sudden growth up to the 
coronal value of 1.5 MK at y = 10 Mm. Then with Eq. (flU)) we obtain the 
corresponding gas pressure and mass density profiles. Both p c (y) (not shown) 
and g c (y) (Fig. [I] right-top panel) experience a sudden fall-off from photosphere 
to coronal values at the transition region. 

2.2. Initial Conditions 

At t = s, we initially perturb the equilibrium impulsively using local- 
ized Gaussian pulses simultaneously in the gas pressure and vertical 
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component of velocity, viz., 



p(x,y,t = 0)=p e (y)+A p f(x,y), V y (x,y,t = 0)=A v f{x,y), (12) 

(x-x Q ) 2 (y-y ) 2 ' 



f(x,y) = exp 



UK. 



(13) 



Here A p and A v are the amplitudes of the perturbations, (xo,yo) is 
their initial position and (w x ,w y ) denotes their widths along the x- 
and y-directions. We set and hold fixed xq = Mm, yo = 0.5 Mm, w x = 0.5 



Mm, w y 



0.5 Mm, A p = 0.02p o (y ) and A v = 0.2 x 10~ i Mm s _1 . These 



magnitudes of A v and A p represent the recently detected flow and 
temperature in solar granulation (Baran, 2011). We separately consider 
three orientations for our unidirectional equilibrium magnetic field lines: 
(a) horizontal, (b) vertical, (c) oblique to the gravitational force. 



3. Numerical results 

Equations ([I}— © are solved numerically using the FLASH code (Lee & Deane, 
2009). This code implements a second-order unsplit Godunov solver with var- 
ious slope limiters and Riemann solvers, as well as adaptive mesh refinement 
(AMR). The main advantage of using AMR technique is to refine a 
numerical grid at steep spatial profiles while keeping a grid coarse at 
the places where fine spatial resolution is not essential. Such AMR 
technique usually introduces interpolation errors at different-size nu- 
merical cells. These errors can result in vertical flow that, albeit 
initially small, can grow with height to an unacceptable magnitude. 
A remedy of this inherent phenomenon is to refine the whole region 
below the transition region, which was adopted in our simulations. We 
used the Roe solver and minmod flux limiter, and set the simulation 
box for the horizontal and vertical equilibrium magnetic field cases in 
the x-direction as 5 < x < 5 Mm. In order to trace plasma structures 
which extend more horizontally in the case of the oblique magnetic 
field we use —2 < x < 8 Mm in this case. Along the y-direction the 
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Figure 2. Numerical blocks with their boundaries (solid lines) and the pulse in velocity 
of Eq. (|12|l (color maps) at t = s for the cases of the vertical and horizontal equi- 
librium magnetic fields. Maximum velocity (red color) corresponds to 0.2 x 1CP 3 
Mm s" 1 . 



numerical box was set as — 1 < y < 19 Mm in all cases. At all boundaries, 
we fix all plasma quantities to their equilibrium values, which lead only to small 
numerical reflections of incident wave signals. Additionally, we increased the 
physical domain size and damp incident waves by adopting a coarse 
numerical grid at the top boundary to avoid influence of reflections 
during the simulation time range. As a result, these reflections do not 
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exert any noticeably effect on the dynamics of the system. In all our studies, 
we use an AMR grid with a minimum / maximum level of refinement set to 
3/7. The initial system of blocks is displayed in Fig. [2j For the vertical and 
horizontal magnetic fields initially, at t = s the whole simulation region is 
covered by 2666 blocks, with a similar number of blocks for the oblique magnetic 
field. As every block consists of 8 x 8 numerical cells, therefore, this number of 
blocks corresponds to 170624 numerical cells. This results in the finest (poorest) 
resolution of Ax = Ay ~ 19.5 x 10~ 3 Mm (Ax = Ay ~ 625 x 10~ 3 Mm) in 
the region below y = 3.7 Mm (above y = 9 Mm) at t = s (see Fig. [5]). The 
refinement strategy is based on controlling numerical errors in mass density, 
which results in excellent resolution of steep spatial profiles and greatly reduces 
numerical diffusion at these locations. The duration of a typical numerical 
run was 1500 s. Although the numerical simulations have been carried out for 
—5 Mm < x < 5 Mm, —1 Mm < y < 19 Mm, only the results in the regions of 
interest are displayed in this work. 

3.1. Vertical equilibrium magnetic field: s = y 

First, we investigate the case of the vertical equilibrium magnetic field, which 
corresponds to § = y in Eq. ([5]). Here y is a unit vector along the y-direction. 
The initial pulse of Eq. (|12[) triggers magnctoacoustic-gravity waves in the solar 
atmosphere. These waves are compressive, and so we can trace their evolution 
in both temperature and velocity. Figure [3] illustrates logarithm of temperature 
profiles (colour maps) and velocity vectors resulting from the perturbations. At 
t = 200 s (top- left panel), the excited perturbations in the vertical compo- 
nent of velocity already penetrated the solar corona arriving to the point 
(x = Mm, y = 5.5 Mm). The perturbation in the vertical component of 
velocity propagates by rj5 Mm in first 200 s of its evolution, which gives 
the speed of propagation as w25 x 10 -3 Mm s _1 . The average sound 
speed in this region is significantly higher. This is also confirmed by 
Fedun et al. (2011b) (cf., their Figs. 4-7) as their waves propagate from 
the transition region up by 2 Mm within about 40 s of the simulation 
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Figure 3. Spatial profiles of logarithm of temperature (color maps) and velocity (ar- 
rows) profiles at t = 200 s (top-left), t = 400 s (top-right), t = 800 s (bottom-left), 
t — 10 3 s (bottom-right), for the case of the vertical equilibrium magnetic field. Temper- 
ature is expressed in units of 1 MK. The arrow below each panel represents the length 
of the velocity vector, expressed in units of 10 x 10~ 3 Mm s _1 . The corresponding 
movie can be found in the file fig3.avi online. 
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time. Prom this comparison we infer that the perturbation we initially 
imposed triggers flows in the solar atmosphere first, and thereafter 
these flows produce magnetoacoustic-gravity waves. At this moment of 
time, the transition region experiences small perturbation and the signal reflected 
from the transition region is represented by arrows around x = Mm, y = 1.5 
Mm. This reflected signal clashes with the slowly up-going wave that results at 
the launching place, i.e., x — Mm, y = 0.5 Mm. As a result of this clash, 
the leading front is generated. Cold chromosphcric plasma lags behind and at 
t = 400 s (Figure [21 top-right panel) it reaches an altitude of y ~ 6.5 Mm. The 
reason for the cold material being lifted is the pressure gradient which works 
against gravity and forces the chromosphcric material to penetrate the solar 
corona (Sterling & Hollweg, 1989). At t = 800 s (Figured bottom- left), the 
transition region experiences well-developed oscillations with the central region 
reaching the level of y ~ 7 Mm at x = Mm, and also the off-side propagating 
oscillations. These oscillations are of lower amplitudes at t = 10 3 s (Figure [21 
bottom- right panel). At this time, the central plasma experiences gravitational 
fall-off while the region of cold plasma widens. The jet, which is associated 
with the shock, lasts for several hundred seconds and its velocity reaches up to 
maximum limit of ~ 60 x 10 -3 Mm s _1 . 




200 400 600 800 
time [s] 

Figure 4. The time signature of V y collected at the point (x — 0, y — 5) Mm for the 
case of the vertical equilibrium magnetic field. 



SOLA: KM-AKS-JAM-R0_Solar_Phys_ver2e.tex; 30 August 2012; 0:11; p. 12 



Magnetoacoustic-Gravity Waves in the Solar Atmosphere 



Figure S] illustrates the vertical component of velocity that is collected in time 
at the detection point (x = 0, y = 5) Mm for the case of Fig. [3] As a result of 
rapid mass density fall off with height, upwardly-propagating waves grow in their 
amplitude and steepen rapidly to form shocks. The arrival of the first shock-front 
at the detection point is clearly seen at t ~ 290 s. The second shock-front reaches 
the detection point at t ~ 540 s, i.e. approximately 250 s later. This secondary 
shock (and subsequent shocks which arrive respectively after about 160 s and 
240 s) results from the nonlinear wake, which lags behind the leading shock 
(Sterling & Hollweg, 1989; Srivastava & Murawski, 2011). These times can be 
compared with the acoustic cut-off period (Roberts, 2006), 

PM= ,/f| y > , (i4) 

which for y = 0.5 Mm attains a value of P ac (y) — 180 s. This value differs 
by 70 s from the time-span between arrivals of neighboring shocks, which is 
250 s. However, the 2D model we discuss here is more complex than 
the ID scenario which is described by the Klein-Gordon equation 
(Roberts, 2006). The wave-period we detected is altered by the interaction 
between up-going waves from the launching place of the initial pulse and the 
reflected one from the transition region signals. As wave reflections result at 
the large temperature gradients, therefore, the chromosphere sustains a cavity 
for these waves which are represented by the oblique stripes located at altitudes 
0.7 Mm < y < 1.4 Mm (Fig. [5]). In the neighborhood of the point x = 0.4 Mm, y = 
0.4 Mm, we observe the formation of vortices (Fig.[5l top panel) which experience 
energy cascade into smaller scales. These vortices are present till the end of our 
simulation runs (Fig. [SJ bottom panel). The first vortex results from the initial 
pulse in V y , which is a characteristic feature of velocity perturbations. This first 
vortex seeds convection in the convectively unstable plasma layers. According 
to the Schwarzschild's instability condition, a medium is convectively unstable 
if the squared buoyancy (or Brunt- Vaisala) frequency, 

UJ > 9 {^-^) = %-{ ' (15) 
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is negative (e.g., Roberts, 2006). As this criterion is satisfied for y < 0.6 Mm, 
therefore, the convection sets in there. Such vortices were also theorized by 
Konkol, Murawski, & Zaqarashvili (2011)in a similar context. 





Figure 5. Total velocity, |V|, (color map) and velocity (arrows) profiles at t = 400 s 
(top) and t — 10 3 s (bottom) for the case of the vertical equilibrium magnetic field. 
The velocity vectors are expressed in units of 1 x 10~ 3 Mm s _1 . The corresponding 
movie can be found in the file fig5.avi online. 

3.2. Horizontal equilibrium magnetic field: § = x 

In this section, we investigate the horizontal equilibrium magnetic field, which 
corresponds to s = x in Eq. (|5|). In this case, the spatial profiles of T drawn 
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Figure 6. Spatial profiles of logarithm of temperature (color maps) and velocity (ar- 
rows) at t = 250 s (top) and t = 10 3 s (bottom) for the case of the horizontal equilibrium 
magnetic field. Temperature is drawn in units of 1 MK. The arrow below each panel 
represents the length of the velocity vector, expressed in units of 2 x 10~ 3 Mm s _1 . 
The corresponding movie can be found in the file fig6.avi online. 

at t = 250 s and t = 10 3 s (Fig. |6|) reveal small amplitude oscillations of the 
transition region without the presence of a jet which was observed for the case 
of the vertical magnetic field. On comparing with Fig. [31 the waves resulting 
from the initial pulses already arrived to the solar corona at t = 250 s (Fig. [6l 
left panel). As the magnetic field lines are horizontal, therefore, these waves 
are essentially fast magnetoacoustic-gravity waves in the region of low plasma 
j3 that occurs for y > 2.6 Mm (Fig. [TJ bottom panel). It is interesting that the 
orientation of magnetic field plays very crucial role in determining the amplitude 
of oscillations of the transition region. 

Similar to the time-signatures of Fig. [U V y collected at the detection point 
(x = 0, y = 5) Mm, reveals the shocks. However, for the horizontal magnetic 
field, we now observe only two shocks (as seen in Fig. [7]). The first shock arrives 
to the detection point at t = 250 s and the second shock reaches this point at 
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Figure 7. The time signature of V y collected at the point (a; = 0, y = 5) Mm for the 
case of the horizontal equilibrium magnetic field. 



t = 400 s. We note that the time-span of 150 s between these events is 100 
s shorter than in the case of the vertical magnetic field in Sect. 13.11 However, 
they are shorter only by 30 s than the acoustic cut-off period, P a c resulting 
from Eq. HU). Note that the subsequent three oscillations, which are of smaller 
amplitude, reveal waveperiods of the order of P ac . 

The chromosphere exhibits essentially similar features to those seen for the 
vertical magnetic field of Sect. 13.11 Both, the vortex motion at x = 0.4 Mm, y = 
0.4 Mm and trapped waves in the solar chromosphere for 0.7 Mm < y < 1.4 Mm, 
are again discernible (Fig. [8]). At t — 400 s (top panel), the vortex is well devel- 
oped. This vortex is similar to that of Fig. |5] (top panel). The flow patterns at the 
line y = 1.4 Mm in Figs. [5] and [5] differ in some details though. In the case of the 
vertical magnetic field, the resultant waves are essentially slow magnctoacoustic- 
gravity waves, while Fig. [8] corresponds to the horizontal background magnetic 
field with much contribution from fast magnetoacoustic-gravity waves. At a later 
time, turbulence results from this vortex and it differs in some details from its 
analog of Fig. [5] (bottom panel) . In both the vertical and horizontal magnetic 
fields, the turbulence originates from the impulsively triggered perturbations 
that are initially launched in the convectively unstable atmospheric layers. 
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Figure 8. Total velocity, |V|, (color map) and velocity (arrows) profiles at t = 400 s 
(top) and t = 10 3 s (bottom) for the case of the horizontal equilibrium magnetic field. 
The velocity vectors are expressed in units of 1 x 10~ 3 Mm s _ . The corresponding 
movie can be found in the file fig8.avi online. 

3.3. Oblique equilibrium magnetic field: s = (x + y) / V% 

In this section, we discuss the case of the oblique magnetic field, for which 
s in Eq. §5§ is at an angle of ir/4 to the horizontal axis. Figure [5] displays 
temperature profiles at t = 400 s (top- left panel), t = 600 s (top-right panel), 
t = 800 s (bottom-left panel) and t = 10 3 s (bottom-right panel). These profiles 
are more complex than those for either the horizontal or the vertical magnetic 
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field of Sect. O and Sect. O At t = 400 s (Figure M top-left), we clearly see 
the rising transition region. At t = 10 3 s (Figure [51 bottom-right), we observe 
surface waves propagating along the transition region. However, for the oblique 
equilibrium magnetic field, we find that time signatures of V y reveal shocks, and 
the waveperiods are about 250 — 300 s (Fig. [TU|) . The jet, which is associated 
with the shock, lasts for several hundred seconds and its velocity reaches up to 
maximum limit of - 30 x 10~ 3 Mm s" 1 . This is a result of complex interaction 
of waves with the background plasma in upper regions of the solar atmosphere. 
This complexity is seen at t = 10 3 s in the velocity profiles of Fig. [Til which 
clearly illustrates the trapped and reflected waves in the solar chromosphere as 
well as well-developed vortices. 

4. Discussion and conclusion 

We have investigated the impulsive excitation of magnctoacoustic-gravity os- 
cillations and have compared and contrasted their resultant propagation under 
different orientations of equilibrium magnetic fields. We performed 2D numerical 
simulation of the velocity and gas pressure pulses, which mimic a solar granule. 
These pulses were initially launched at the top of solar photosphere, in a stratified 
solar atmosphere utilizing the VAL-C temperature profile. We find that in the 
cases the background magnetic field possesses a non-zero vertical component 
the amplitude of the upwardly-propagating perturbations rapidly grows with 
height due to the rapid decrease in the equilibrium mass density. Therefore, 
the perturbation quickly steepens into shocks in the upper regions of the solar 
chromosphere that launches the cool material behind it. 

We find that the excitation of magnetoacoustic-gravity waves in a non-ho- 
rizontal equilibrium magnetic field configuration is channeled upwards and is 
able to cause large amplitude transition region oscillations. Meanwhile, the hor- 
izontal equilibrium magnetic field configuration results in comparatively smaller 
amplitude transition region oscillations, due to fast magnetoacoustic-gravity 
waves which spread their energy across magnetic field lines in contrast to slow 
magnctoacoustic-gravity waves which in a strongly magnetized plasma region 
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Figure 9. Spatial profiles of logarithm of temperature (color maps) and velocity (ar- 
rows) profile at t = 400 s (top- left), t = 600 s (top-right), t = 800 s (bottom- left), 
t = 10 3 s (bottom-right) for the case of the oblique equilibrium magnetic field. Tem- 
perature is drawn in units of 1 MK. The arrow below each panel represents the length 
of the velocity vector, expressed in units of 4 x 10~ 3 Mm s . The corresponding 
movie can be found in the file fig9.avi online. 
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Figure 10. The time signature of V y collected at the points (a; = 4.5, y = 5) Mm (top) 
and (x — 3.84, y = 4.98) Mm (bottom) for the case of the oblique equilibrium magnetic 
field. 




■ 



Figure 11. Total velocity, |V|, (color map) and velocity (arrows) profile at t = 10 3 s for 
the case of the oblique equilibrium magnetic field. The velocity vectors are expressed 
in units of 1 x 10~ 3 Mm s _1 . The corresponding movie can be found in the file 
figll.avi online. 
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are guided along magnetic field lines. This indicates that the wave energy is 
transfered into the upper atmosphere in larger amounts where the magnetic 
field is more vertical. However, the complexity and therefore the evolution of a 
horizontal equilibrium magnetic field allows the reflection and trapping of the 
waves in the lower solar atmosphere and can influence the localized dynamics 
and heating of the atmosphere. The reflection of the waves in the lower solar 
atmosphere and its trapping may also result in the presence of chromospheric 
cavities. 

By analyzing the time-signature of V y collected at a fixed spatial point, we 
also find that the period of oscillation is longer when comparing the vertical 
magnetic field case (here, the characteristic period of oscillation was « 250 — 300 
s) with the horizontal equilibrium magnetic system (in which the characteristic 
period of oscillation was « 150 — 200 s). For the vertical magnetic field system 
in low plasma /3 regions, the magnetoacoustic-gravity waves are well-described 
as slow magnetoacoustic-gravity waves (propagating along the magnetic field 
lines at approximately the sound speed, c s (y), see Eq. [7]). For the horizon- 
tal magnetic field system in the strongly magnetized plasma, the oscillations 
transverse to the magnetic field are well-described as fast magnetoacoustic- 
gravity waves (propagating across the magnetic field lines at the fast speed, 
Cj (y) = c\(y) + Cg(y), see Eqs. |5]and[7]). As for y = 0.5 Mm we have ca < c s , 
therefore the cut-off frequencies of the fast and slow magnetoacoustic-gravity 
waves are very close to each other, which results in very close values of the cut- 
off waveperiods. Therefore, we could expect that the detected waveperiods for 
the vertical and horizontal magnetic fields exhibit similar values. However, the 
cut-off waveperiods are derived on the basis of the linear theory which is valid 
for small amplitude oscillations only. As such small oscillations are present in the 
case of the horizontal magnetic field, therefore, the obtained numerical data lies 
close to the analytical prediction. In the case of the vertical magnetic field, the 
upwardly propagating waves interact with the waves which become reflected 
from the transition region. As a result of that larger amplitude oscillations 
originate, which significantly alter the background plasma. The waves, which 
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propagate through this strongly modified medium exhibit modified velocities 
and they get reflected from the transition region that is locally largely curved. 
As a consequence of that waveperiods within the range 250 — 300 s result, which 
differ from P ac . 

It is known that radiation is an effective mechanism of wave damp- 
ing in the low photosphere ((Mihalas &: Toomre, 1982). It might not 
radically alter the system's behavior, but radiative damping is at least 
likely to reduce the amplitude of the waves reaching the transition 
region, leading to shorter jets than what we see in our models, as 
well as lower amplitudes of the coronal shocks. In this paper, we do 
not invoke the radiative cooling nor thermal conduction in our model 
atmosphere as we aim to model the small-scale atmospheric regions 
above the solar photosphere where such effects are not believed to be 
dominant. However, we intend to include these effects in our future 
studies. 

It is noteworthy that the solar atmosphere is structured by convec- 
tive overshoot which is absent in the model we devised. Instead, we 
isolated a single granule-like perturbation by aiming to mimic a mag- 
nitude of flow and plasma temperature that are associated with the 
solar granulation (Baran 2011), and considered the complex scenario 
which results by this simple model to understand the physics of wave 
phenomena above such magnetic structures in the solar atmosphere. 
We intend to develop more advanced models in our future studies. 

In conclusion, our numerical simulations clearly demonstrate that small am- 
plitude initial pulses in vertical velocity and gas pressure are able to trigger a 
plethora of dynamic phenomena in the upper regions of the solar atmosphere 
with waveperiods within the range of 150 — 300 s, a value which depends on 
orientation of the background magnetic field. However, it should be noted that 
the performed 2D simulations are idealized in the sense that they do not include 
radiative transfer and thermal conduction along field lines. The magnetic field 
configuration and the equilibrium stratification are simple and we modeled a 
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single granule only. These limitations require additional studies which we intend 
to carry on in near future. 
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